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Aims. We studied the multi-frequency predictions of various annihilating Dark Matter (DM) scenarios in order to explore 
the possibility to interpret the still unveiled origin of non-thermal phenomena in galaxy clusters. 

Methods. We consider three different DM models with light (9 GeV), intermediate (60 GeV) and high (500 GeV) 
neutralino mass and we study their physical effects in the atmosphere of the Coma cluster. The secondary particles 
• created in the neutralino annihilation processes produce a multi-frequency Spectral Energy Distribution (SED) of non- 

thermal radiation, as well as heating of the intracluster gas, that are tested against the observations available for 
the Coma cluster, from radio to gamma-rays. The various DM produced SEDs are normalized by the condition to fit 
, the Coma radio halo spectrum, thus obtaining best-fit values of the annihilation cross-section aV and of the central 

magnetic field Bo- 

Results. We find that it is not possible to interpret all the non-thermal phenomena observed in galaxy clusters in terms 
of DM annihilation. The light mass DM model with 9 GeV mass produces too small power at all other frequencies, 
while the high mass DM model with 500 GeV produces a large excess power at all other frequencies. The intermediate 
mass DM model with 60 GeV and composition is marginally consistent with the HXR and gamma-ray observations 
^ ' but barely fails to reproduce the EUV and soft X-ray observations. The intermediate mass DM model with 60 GeV 

\^ ' and bb composition is, on the other hand, always below the observed fluxes. We note that the radio halo spectrum of 

00 . Coma is well fitted only in the bb or light and intermediate mass DM models. We also find that the heating produced 

^SJ ■ by the DM annihilation in the center of the Coma cluster is always larger than the intracluster gas cooling rate for an 

NEW DM density profile and it is substantially smaller than the cooling rate only for a cored DM density profile in 
light mass DM model with 9 GeV. 

Conclusions. We conclude that the possibility of interpreting the origin of non-thermal phenomena in galaxy clusters 
with DM annihilation scenarios requires a low neutralino mass and a cored DM density profile. If we then consider the 
multimessenger constraints to the neutralino annihilation cross-section, it turns out that such scenario would also be 
excluded unless we introduce a substantial boost factor due to the presence of DM substructures. 
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1. Introduction al 1999), even if sometimes the observations are controver- 

. , . , , . , sial (e.g. Rossetti & Molendi 2004). 
Observational evidence tor the existence oi dmuse compo- 
nents of non-thermal origin in clusters of galaxies and their While the diffuse radio sources in galaxy clusters are 
wider implications is a topic of sustained debate among universally interpreted as produced by synchrotron emis- 
researchers. In particular, the presence of giant radio ha- sion, the Extreme Ultraviolet (EUV) and Hard X-ray 
los and relics (see Feretti & Giovannini 2008 for a review) (HXR) excesses have been interpreted as the combined ef- 
indicates that diffuse distributions of relativistic electrons feet of an underlying inverse-Compton scattering (ICS) of 
and magnetic fields must exist in many clusters of galaxies, the Cosmic Microwave Background (CMB) against cluster 
Moreover, some clusters show also the presence of 'unusual' relativistic electrons that complement the radio halo popu- 
diffuse components, i.e. the excess soft and hard X-ray emis- lation (e.g., Hwang 1997, Ensslin & Biermann 1998, Sarazin 
sion above the level expected from the hot virialized cluster & Lieu 1998, Blasi & Colafrancesco 1999, Colafrancesco 
(e.g.. Lieu et al 1996, Fusco-Femiano ct al 1999, Kaastra et et al. 2005, Marchegiani, Perola & Colafrancesco 2007, 

Colafrancesco & Marchegiani 2009, Buote 2001), or ther- 
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'missing baryons' at sub-virial temperatures (see, e.g., the 
review of Durret et al 2008). Focusing upon the former 
interpretation, the role and ramifications of an underlying 
power law spectrum across the entire X-ray frequency band 
can be far reaching, with issiics concerning cluster merging 
and central cooling, and also mass and spectroscopic conse- 
quences being addressed by Million & Allen (2009), Lagana 
et al (2010), and Prokhorov (2009). In the even broader 
context, the observed cosmic ®Li abundance could be evi- 
dence for excess entropy in clusters (Nath ct al. 2006), along 
with the surprisingly low level of Sunyaev-Zel'dovich effect 
in the WMAP and SPT sample of many clusters (Diego 
and Partridge 2010, Lieu et al. 2006, Komatsu et al. 2010, 
Luckcr ct al. 2009), and the effects of the magnetic field 
on the intracluster gas distribution (e.g. Colafrancesco & 
Giordano 2007) as possible indication of significant intra- 
cluster non-thermal pressure. 

Historically however, there have been numerous sugges- 
tions of a general, non-thermal intracluster environment, 
early examples being Jaffe (1977) and Rephaeli (1979) who 
considered Alfven wave excitation. In particular, a model of 
the acceleration via Alfven waves driven by major cluster 
mergers is given in Brunetti et al. (2004), including the idea 
that relativistic jets carry a significant portion of the total 
energy output of radio galaxies, causing X-ray emission up 
to Mpc distance scales (Ghisellini & Celotti 2001). Celotti 
ct al. (2001) provided an analysis of such a scenario for the 
case of PKS 0637-752. 

A relevant, related observation (see Giovannini et al. 1993) 
maps the spectral index of cluster radio halos, and specif- 
ically of the Coma radio halo from 327 MHz to 1.4 GHz. 
Here, we note the need for non-thermal emission by elec- 
trons with energy just above the energy range required by 
the non-thermal electron population supposed to produce 
soft and hard X-rays by ICS (see Lieu et al. 1999). 

An important issue connected with the presence of 
diffuse emission of non-thermal origin in galaxy clusters is 
that relativistic electrons with E ~ GeV (as those which 
produce the radio halos) lose their energy in relatively 
short time by effect of synchrotron and ICS with the 
CMB photons interactions, and they should reach a short 
distance in comparison to the dimension of the halos (see 
discussion in Blasi & Colafrancesco 1999). Therefore, sev- 
eral solutions have been proposed to solve this problem, as 
the re-acceleration in situ of the electrons by Alfven wave 
excitation driven by major cluster mergers and induced 
turbulences (e.g., Jaffe 1977, Rephaeli 1979, Brunetti et al. 
2004) or the production of secondary electrons in the in- 
teractions between non-thermal and thermal protons (e.g., 
Dennison 1980, Blasi & Colafrancesco 1999, Marchegiani 
et al. 2007, Colafrancesco & Marchegiani 2008). 

The separate physical mechanism to be invoked is 
ncutralino Dark Matter (DM) annihilation as a cluster's 
non-thermal reservoir (sec Colafrancesco & Mele 2001, 
Colafranceso, Profumo, & Ullio 2006) to compensate for 
losses of the particles as they are distributed throughout 
the cluster medium during the long process of particle dif- 
fusion (see, e.g., Volk 1996, Colafrancesco & Blasi 1998). 

Today, several astrophysical evidences (e.g., gravita- 
tional lensing, galaxy rotation curves, galaxy clusters 
masses) indicate that most of the matter content of the 
universe is in form of DM, whose natiirc is still elusive. 
Several candidates have been proposed as fundamental DM 



constituents, ranging from axions to light, MeV DM, from 
Kaluza-Klein (KK) particles, branons, primordial Black 
Holes (PBH), mirror matter to supersymmetric Weakly 
Interacting Massive Particles (WIMPs; see, e.g., Baltz 2004, 
Bertone et al. 2005, and Bergstrom 2000 for reviews). In 
this paper we will assume that the main DM constituent 
is the lightest neutralino of the minimal supersymmet- 
ric extension of the Standard Model (MSSM). This as- 
sumption can be tested by observing the radiation emit- 
ted by the secondary products of the neutralino annihi- 
lation processes, whose features have been studied (e.g., 
Colafrancesco & Mele 2001, Colafrancesco et al. 2006). As 
pointed by Colafrancesco et al. (2006), the relevant physi- 
cal properties which determine the features of the emitted 
radiation are the composition of the neutralino, its mass, 
and the value of the annihilation cross section. 

In this paper we explore the consequences of a neu- 
tralino DM annihilation scenario on the non-thermal parti- 
cle content of galaxy clusters over a multi-frequency range 
(from radio to gamma-ray and TeV energies), and in partic- 
ular on one of the best studied clusters so far, like the Coma 
cluster of galaxies, for which an extensive multi- frequency 
coverage has been accumulated. 

To this aim, we focus on three neutralino DM models with 
mass of 9, 60 and 500 GeV. 

Neutralinos with 60 and 500 GeV mass are representa- 
tive of, respectively, intermediate and high mass supersym- 
metric models which are consistent with accelerator and 
multi-messenger constraints, and that are also interesting 
to be probed by direct detection experiment (like e.g. the 
recent CDMS results, see Ahmed et al. 2009). In this con- 
text, the 60 GeV neutralino mass is close to the lower limit 
of the combination of the CDMS II experiment with the 
allowed parameter space calculated from various Minimal 
Supersymmetric Models (e.g., Roszkowski et al. 2007, Ellis 
et al. 2005), and it is interesting for being possibly tested 
by a multi-frequency search limited to the gamma-ray ener- 
gies probed by Fermi. On the other side, the 500 GeV mass 
is interesting for being tested by a multi-frequency search 
extended even in Cherenkov experiments realm. 

The light mass neutralino model with = 9 GeV 
has been recently worked out as a consistent solution of 
the CoGeNT and DAMA observations (Fitzpatrick et al. 
2010). Such a low mass neutralino is allowed in supersym- 
metric scenarios where the hypothesis of gaugino masses 
universality at the unification scale is released (Bottino et 
al. 2003). 

To explore the astrophysical consequences of such DM 
models, we consider pure annihilation channels in order to 
highlight their observational features. We focus, specifically, 
on bb and annihilations for 9, 60 and 500 GeV mass 
neutralino models, and we also consider annihilations into 
for the 500 GeV mass ncutralino model. We leave the 
neutralino annihilation cross-section aV as a free parameter 
to be constrained, in our analysis, by the astrophysical data. 
In this study, we also consider the constraints on the DM 
model and density profile coming from the heating rate of 
the thermal gas at the centre of the cluster. 

The outline of the paper is the following. We describe in 
Sect. 2 the properties of the various neutralino DM models 
and the various DM density profiles that we consider for the 
Coma cluster. We derive the astrophysical consequences of 
the neutralino DM annihilation in Coma in Sect. 3, and we 
focus on the limits set on the expected DM-induced spec- 
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tral energy distribution (SED) at niultifrequency, as well 
as on the additional physical effects induced by DM anni- 
hilation. The discussion of our results and the conclusion 
of our analysis are presented in Sect. 4. 

Throughout the paper, we use a flat, vacuum- 
dominated cosmological model with fi^ = 0.3, ilA = 0.7 
and Ho = 70 km s"^ Mpc'^. 



profile a ( kpc) — 

NFW 1 287 1.96 X 10* 

cored 142 1.66 x 10^^ 

NFW^-''"'''^ 1 450 6.75 x 10^ 

Table 1. Density profiles for the Coma cluster, pc = 

1.36 x 10^MQkpc~ is the critical density for closure of 
the universe. 



2. The Dark Matter model 



A 9 or 60 GeV supersymmetric particle is expected to anni- 
hilate in quarks and leptons. Direct annihilation into light 
fermions is helicity suppressed, so that the main remain- 
ing channels of annihilations are bb and t~^t~. Electrons 
and positrons will be, then, mainly produced through sec- 
ondary processes of hadronization and/or decay. Direct an- 
nihilation into 7S happens at 1-loop level and is therefore 
highly suppressed too. The bulk of photons arises from the 
^0 _^ g decay. 

In addition to the aforementioned channels, a 500 GeV 
supersymmetric particle can also annihilate into gauge 
bosons, heavy quarks and Higgs particles. The yield of pho- 
tons produced in these cases, as well as their energy spec- 
trum, is very much similar to each other, so that we choose 
as a representative case the annihilation into W^W~. 

In this framework, we have computed the 
positron/electron spectra with DarkSUSY (Gondolo 
et al. 2004) and the yield of gamma-ray photons according 
to Fornengo, Fieri & Scopel (2004). 

2.1. The Dark Matter profile of the Coma cluster 

Recent determinations of the mass and density profile of 
Coma include methods based on weak lensing (Kubo et al. 
2007) as well as on the study of the velocity moments of 
early-type galaxies in the Coma cluster (Lokas & Mamon 
2003) or on the analysis of infall patterns around Coma 
(Geller et al. 1999). Numerical simulations also give a hint 
on the DM density profile inside clusters (Bullock et al. 
2001). The results of astronomical measurements and nu- 
merical experiments lead to a determination of the mass 
profile which is consistent among the various methods, 
within the errors and the spatial resolution of each ap- 
proach. In this paper we use the results of Lokas & Mamon 
(2003), since they provide a fit to the data with both a 
spiky NFW profile and a cored profile. 

The total inferred dark matter mass of Coma is M^r — 
1.2 • 10^^ Mq within a dark matter virial radius of 2.7 Mpc. 
The profile is parametrized as: 



p{r) 



1 



(1) 



Scale parameters and inner slope can be found in Table 
[TJ Also shown are the parameters for an NFW profile as 
inferred by N-body simulations (Bullock et al. 2001). 

We also compute the effective boost factor seen as a line 
of sight (los) integral due to the presence of subhalos orbit- 
ing in the gravitational potential of the Coma cluster. We 
follow the methodology and results of Fieri et al. (2009) in 
the framework suggested by the Via Lactea II numerical 
simulation (Diemand et al. 2008). 

We make therefore the following assumptions: i) the slope 



^ 10-= 



smooth I 
su b halos 



smooth 



10-1 
0.01 



1 

O(dcg) 



Fig. 1. The los integral of the DM density squared is shown as 
a function of the projected angle from the cluster center for a 
resolution of 0.1 degree (corresponding to a solid angle of 10~^ 
sr). The red dotted curve is for a smooth DM density profile 
while the blue sohd curve is for a smooth DM density profile 
with the addition of subhalos (see text for details). 



of the subhalo mass function is -2, which implies, for a DM 
halo such as Coma, that about 64% of its mass is virial- 
ized in subhalos; ii) the mass of subhalos is in the range 
1O~^M0 — 10~'^ Myir] in) the radial dependence of the sub- 
halo concentration parameter; iv) an NFW profile for each 
of the subhalos. 

We then compute the los integral, with and without subha- 
los, of the DM density squared in order to infer the boost 
factor B defined as the integral along the line of sight and on 
the solid angle of the neutralino density squared in the case 
in which we consider a smooth DM density profile or the 
sum of smooth DM density plus DM subhalos (see, e.g.. 
Fieri et al. 2009). Fitting formulas for the concentration 
parameter as well as explicit notation for the line of sight 
integral including a population of subhalos can be found in 
Fieri et al. (2009). 

FigH] shows the result of our computation in the case of 
an angular resolution AJ7 = 10^^ sr (corresponding to 0.1 
degree). Once integrating over the whole cluster volume, 
we obtain a boost factor B — 35. We do not expect signif- 
icant deviation from this value when changing the subhalo 
profile. 
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In the following calculations we will first consider the 
case of a smooth DM density profile and then the effects of 
the boost factor. 

3. Astrophysical consequences 

We discuss here various astrophysical consequences of the 
DM models we consider in this paper: e.m. signals (all of 
non-thermal origin) and heating of the intracluster gas in- 
duced by secondary particles produced in the DM annihi- 
lation processes. 

The complete description of the emission features in- 
duced by DM annihilation is done by using the diffusion 
equation (i.e. neglecting convection and re-acceleration ef- 
fects) given in Colafrancesco, Profumo & Ullio (2006, see 
their eq. 35) which takes into account, consistently, the dif- 
fusion and energy loss properties of the secondary particles 
produced in the neutralino DM annihilation process. 

The energy loss term is the sum of effects due to 
Inverse Compton, synchrotron radiation. Coulomb losses 
and Bremsstrahlung (see eq. A. 16 in Colafrancesco et al. 
2006 for details). For electron energies ^ GeV the Inverse 
Compton and synchrotron terms dominate the losses while 
for energies below ^ 140 MeV, Coulomb losses dominate 
(see Fig. A. 3 in Colafrancesco et al. 2006). For a more de- 
tailed description of the energy losses in the Coma clus- 
ter, we use throughout the paper a B-field radial profile 
B{r) oc nthif), where nth is the intracluster gas number 
density, unless otherwise specified. 

The derivation of the full solution of the diffusion equa- 
tion and the effects of diffusion and energy losses are given 
in Colafrancesco et al. (2006, see Eq. A.l) and we refer 
to this paper for technical details. We consider here the 
same analytical approach and we present the results of the 
multi-frequency emission produced by the DM models un- 
der study for the Coma cluster. We also compute the DM- 
induced signals in the limit in which electrons and positrons 
lose energy on a timescale much shorter than the timescale 
for spatial diffusion, i.e. the regime which applies to the 
case of galaxy clusters (see Colafrancesco et al. 2006 for a 
discussion) . 

The part of the SED from which the strongest spectral 
constraints can be set is the low-frequency radio band from 
~ 30 MHz to - 5 GHz. 

FiglH shows the spectra of the diffuse radio emission of 
Coma calculated in the various DM models that we con- 
sider here. The best fit values of the central value of the 
magnetic field Bq and aV have been obtained from the 
fit to the Coma radio halo spectrum and normalizing the 
curves at the data point at 608 MHz. 

We notice that the radio-halo spectrum could discriminate 
between the different neutralino compositions, preferring a 
bb composition due to its better-fitting intrinsic spectral 
curvature (see Figl2). 

The fit is influenced by the value of the B-fleld. 
Changing the intensity of the B-fleld implies to change the 
energies of the electrons that emit by synchrotron in this 
frequency range, and this implies, in turn, that the result- 
ing radio halo spectrum has a change in the spectral shape 
according to the value of Bq- These effects are illustrated 
in Figl21 For = 9 GeV, the bb model recovers quite well 
the curvature of the observed spectrum for Bq — 20 ^G, 
while the model reproduces the observed spectrum for 
a value of Bq = 2.2 ^G. For = 60 GeV, the bb model 



recovers quite well the curvature of the observed spectrum 
for Bq — 0.5 fj,G, while the model reproduces the ob- 
served spectrum only for the frequency range v > 100 MHz 
for a value of Bq = 0.1 fiG. For = 500 GeV, the bb 
model reproduces the curvature of the observed spectrum 
for Bo = 0.01 //G, while the and models require 
values of the central magnetic fleld that are much lower 
{Bq = 0.002 — 0.001 /iC) and can reproduce the spectral 
shape of the radio halo only for > 400 MHz. The cor- 
responding best fit values of the annihilation cross section 
aV are given in the caption of Figl^l 

FiglS] shows the radio halo spectrum calculated in the 
various DM models that we consider here where we assumed 
the value of the central magnetic field derived from the most 
recent analysis of Faraday Rotation measures, Bq — 4.7 fiG 
with a radial profile given by B{r) oc nthi'f'Y''^ (Bonafede et 
al. 2010). In this figure we have fixed the value of aV such 
that the radio halo spectrum fits the observed one at 608 
MHz. For this value of the central B-field, the DM model 
with — 9 GeV and composition and the model with 
= 60 GeV and bb composition can (marginally) repro- 
duce the observed shape of the spectrum, while the model 
with 9 GeV and bb composition is too steep and all other 
DM models provide radio halo spectra much flatter (harder) 
than the measured one. 

The overall spectral energy distributions (SEDs) of the 
DM models derived from the fits to the Coma radio halo 
spectrum populate differently the high-frequency regions of 
the e.m. spectrum. 

The SEDs of a 9 GeV, of a 60 GeV and of a 500 GeV 
neutralino models annihilating in Coma are shown in Fig. 

SI 

The 9 and 60 GeV mass DM models that best 
fit the radio-halo spectrum of Coma produce ICS and 
bremsstrahlung radiation and 7r° — > 77 radiation consis- 
tent with the existing limits at multi-frequency. Under the 
previous conditions, these signals cannot, in particular, be 
the nature of the soft and hard X-ray emission excesses de- 
tected in this cluster. Only the 60 GeV neutralino model 
can fit both the radio halo spectrum and the hard X-ray 
data, but not the Extreme UV Explorer (EUVE) data. 

The DM model with 500 GeV mass can reproduce the 
Coma radio halo spectrum only with very low values of the 
central magnetic field, Bq ~ 0.001 — 0.1 /iC. This fact has 
the consequence of requiring quite large values for the anni- 
hilation cross section aV ~ 10~^^ — 10~^^cm^/s, and hence 
for the secondary electron density produced by the neu- 
tralino annihilation. Therefore, we obtain that ICS emis- 
sion expected from these models is very large and exceeds 
both the EUVE and BeppoSAX data as well as the EGRET 
and Fermi upper limits on Coma (the HESS upper limit on 
Coma at very high energies cannot strongly constrain the 
SED produced in the DM models that we consider here). 
We conclude, therefore, that neutralino DM models with 
mass of order of 500 GeV must be excluded by our analy- 
sis. 

The shape of the SED does not change substantially 
when a cored DM profile is considered. The only effect is 
that the cored DM profile requires slightly larger values on 
aV to fit the Coma radio halo data. 

We also evaluate the intracluster gas heating induced 
by DM annihilation (see Figl5]) and we find that an NEW 
DM density profile produces an excess heating in the inner 
^ 3 kpc and ~ 30 kpc regions of Coma for 9 and 60 GeV, 
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respectively. This effect provides, hence, a strong constraint 
to the annihilation cross-section for these DM models. In 
order to have the DM-induced heating rate lower than the 
bremsstrahlung cooling rate at the cluster center the anni- 
hilation cross section worked out here must be reduced by 
a factor - lO'^ - lO^^ 

The heating rate for the case of a cored DM profile 
with a mass of 500 GeV is larger than the bremsstrahlung 
cooling rate of Coma for every DM annihilation channel. 
Therefore, also this result excludes the 500 GeV mass neu- 
tralino model. 

DM models with a mass of 60 GeV and a cored profile show, 
instead, a heating rate smaller than the bremsstrahlung 
cooling rate for the 66 composition, and are therefore consis- 
tent with both multi-frequency SED data and with heating 
rate constraints, while the model is not consistent be- 
cause it has a heating rate greater than the cooling rate. For 
a neutralino mass of 9 GeV and cored DM profile, both 66 
and models are consistent with heating rate constraints. 



4. Discussion and conclusions 

Any DM interpretation of astrophysical data has intrin- 
sic multi-frequency correlations and should, therefore, be 
tested against a series of independent observational con- 
straints. 

The full SED induced by DM annihilation in Coma in- 
dicates that the non-thermal e.m. signals induced by a 9 or 
60 GeV neutralino can be consistent with the shape of the 
radio-halo spectrum and produce ICS, bremsstrahlung ra- 
diation and TT*^ — >■ 77 radiation consistent with the existing 
limits at multi-frequency. However, these signals cannot, in 
particular, be the nature of the soft and hard X-ray emis- 
sion excesses detected in the Coma cluster. An exception to 
this point is given by the model with mass 60 GeV that 
is able to fit the hard X-ray emission observed in Coma. 
However, in the framework of this model, the EUV excess 
in Coma should then be interpreted as due to thermal emis- 
sion from warm baryons at sub-virial temperature. 

Focusing on the radio frequency range,, we have shown 
that a DM particle annihilating into 66 is best fitting the 
spectral shape of the radio halo of the Coma cluster, quite 
independently of the neutralino mass. However, the radio 
halo spectrum fit requires that i) the value of the neutralino 
annihilation cross section aV strongly increases with in- 
creasing neutralino mass, and that ii) the best fit value of 
the central magnetic field in Coma decreases with increas- 
ing neutralino mass. 

In the 9 GeV case the value of the central magnetic field 
required in the 66 case {Bq ^ 20 ^G) is, however, too 
high compared to the results of Faraday Rotation measures 
(CariUi & Taylor 2002, Bonafede et al. 2010), while in the 

case an acceptable value of Bq ^ 2.2 /iG is required. 
The best fit values of Bq decrease then to still marginally 
acceptable values ^ 0.1 — 0.5 jiG for the 60 GeV case (still 
consistent with the available observational limits, see Carilli 
& Taylor 2002), while they take extreme values as low as 
0.001 — 0.01 jiG for the highest neutralino mass (i.e. 500 
GeV) we consider in our study. 

We have also studied the variation of aV for increasing 
values of Bq that could still provide reasonable fit to the 
Coma radio halo spectrum. FigjB] shows that the curves in 
the aV — Bq plane that are locus of the best fit to the 



Coma radio halo fiatten to approximately the same value 
of (jV for high values of Bo, i.e. when the energy losses are 
dominated by synchrotron losses. 

The shaded bands in Fig. [6] encompass the range of up- 
per limits on the annihilation cross-section derived from 
multi-messenger constraints. Such constraints, summarised 
in Table [2j include positrons, antiprotons, radio emission 
and 7-rays from the Galactic Centre and the optical depth 
of CMB photons. The uncertainty related to the propaga- 
tion of antimatter in our Galaxy is taken into account by 
using minimal, mean and maximal propagation setups as in 
Donato et al. (2004). Galactic DM distribution is modelled 
according to the smooth and clumpy distributions found in 
the latest high-resolution iV-body simulations, namely Via 
Lactea II (Diemand et al. 2008) and Aquarius (Springel et 
al. 2008). As a conservative case, also a cored isothermal 
profile with no substructure is considered. For further de- 
tails we refer to Pato et al. (2009), Fieri et al. (2009) and 
Catena et al. (2009). We stress that these multi- messenger 
bounds are quite robust and not easily avoidable, and there- 
fore they set further constraints to the seven benchmark 
DM models studied in the present work (see Fig. |6]). Other 
upper limits on the annihilation cross-section may be de- 
rived e.g. from the observation of dwarf galaxies with the 
Fermi satellite (Farmer 2009; see Fig.3 of Abdo et al. 2010), 
which exclude values greater than ~ 10^^*^ cm^s~^ for a 
neutralino mass of 60 GeV. However, taking into account 
the uncertainty in the density profile as well as the de- 
pendence of the Fermi-LAT sensitivity from the assumed 
source spectrum (which is rather steep for the case of the 
neutralino DM annihilation), such a limit could be released 
toward slightly higher values. 

The fiattening of the curves in the aV — Bq plane for high 
values of Bq means that it is not possible to have values of 
aV that match the allowed region of annihilation cross sec- 
tions suggested by multi-messenger studies of DM annihila- 
tion simply by increasing the value of the central magnetic 
field in galaxy clusters. 

The lower set of curves in Fig|6] shows the annihilation 
cross sections required to fit the radio halo spectrum if DM 
substructures in Coma are considered (according to our dis- 
cussion in Sect. 2). The effect of substructures is to increase 
the surface brightness intensity of the radiation emitted by 
the electrons, especially at large radii. As shown in FigEl 
the inconsistency between fitting the radio halo and the 
multi-messenger analysis is significantly improved, partic- 
ularly for Bq ^ 1/iG. However, none of the DM models 
under study is safely below the annihilation cross section 
upper limits obtained from the multi-messenger constraint 
analysis whose results are summarized in Tab HI 

We have also considered the predictions of the various 
DM models by using a B-field of Coma that has the inten- 
sity and radial profile suggested by recent Faraday Rotation 
measures (Bonafede et al. 2010, see Fig. [3]). In this case, the 
DM model with 60 GeV and 66 composition and the model 
with 9 GeV and composition would produce a radio 
halo spectrum with a shape that is still consistent with the 
observed one. The value of the annihilation cross section 
required in the 60 GeV, 66 case is aV = 1.0 x lO"^** cm"^ 
s~^, while in the 9 GeV, t"^ case the cross section required 
is aV = 1.4 X 10-25 cm^ s'^ 

However, in order to have a DM-induced heating rate lower 
than the bremsstrahlung cooling rate of the intracluster gas 
at the cluster center, the annihilation cross sections worked 
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upper limit on the annihilation cross-section (10 cm'^s 



DM benchmark 


VL2 


An 






= 9 GeV 


XX 


r+T- 


(0.1,0.1,0.1) (rad,rad,rad) 


(7.8,7.8,7.8) (opt,opt,opt) 


(7.8,7.8,7.8) (opt,opt,opt) 


Mx 


= 9 GeV 


XX 


bb 


(0.06,0.06,0.06) (rad,rad,rad) 


(4.3,0.7,0.2) (opt,anp,anp) 


(4.3,2.6,0.8) (opt,anp,anp) 


Mx 


= 60 GeV 


XX 


r+T- 


(2.2, 2.2, 2.2) (rad,rad,rad) 


(5.3,3.6,2.1) (pos,pos,pos) 


(12,8.5,8.0) (pos,pos,pos) 


Mx 


= 60 GeV 


XX 


bb 


(0.6, 0.6, 0.6) (rad,rad,rad) 


(8.5, 0.8, 0.3) (anp,anp,anp) 


(19,3.8,1.7) (anp,anp,anp) 


Mx 


= 500 GeV 


XX 


r+r- 


(106, 101, 73) (rad,pos,pos) 


(125,64,34) (GCg,pos,pos) 


(349, 184, 170) (pos,pos,pos) 


Mx 


= 500 GeV 


XX -5- 


w+w 


(20, 20, 7.0) (rad,rad,anp) 


(74,8.6,3.5) (pos,anp,anp) 


(162,46,20) (pos,anp,anp) 


Mx 


= 500 GeV 


XX 


bb 


(12, 12, 7.4) (rad,rad,anp) 


(54,9.2,3.7) (pos,anp,anp) 


(118,49,21) (pos,anp,anp) 



Table 2. The upper limits on the annihilation cross-sections in units of 10"^® cm^s^^ from several measurements for the seven 
benchmark DM models. Each column shows the limits derived using different galactic DM distributions: Via Lactea II (VL2), 
Aquarius (Aq) and a cored isothermal profile (iso). Each cell features three values corresponding to the minimal, mean and 
maximal propagation setups respectively. Also shown is the channel that produces the quoted upper limit — "rad" , "pos" , "anp" , 
"GCg" and "opt" stand for radio, positrons, antiprotons. Galactic Centre 7-rays and CMB optical depth, respectively. The extreme 
values for each benchmark are indicated in bold and delimit the shaded regions of Fig[6l 



out here must be reduced by a factor ~ 10^ — 10^^ for 
from 9 to 500 GeV, if we assume an NFW density profile. 

This problem is, however, not present for a neutralino 
mass of 9 GeV (with both hb and compositions) and 60 
GeV and composition bb if DM density has a cored profile 
(see FiglS]), while it still holds for the 60 GeV mass with 
composition and for all neutralino compositions in the 
500 GeV mass model. 

The high-mass DM model with 500 GeV, whose 

SED is normalized in order to fit the Coma radio-halo data 
with a very low magnetic field value (i?o ~ 0.01 /xG for the 
bb composition), is inconsistent with the multi- frequency 
SED because it exceeds by large factors the EUV, soft X- 
ray, hard X-ray and gamma-ray (in the energy range of 
Fermi) upper limits (see FigH]). For such high- mass DM 
model, also the and cases barely fail to fit the Coma 
radio-halo spectrum at low frequency below ^ 400 MHz 
(see FigE]). 

We note finally that the high-energy HESS gamma-ray 
limit is not relevant in any of the DM models considered in 
this paper, because in all the cases the threshold energy of 
the HESS upper limits is larger than the neutralino masses 
here considered. 

As a final remark, we stress that the DM candidates 
with = 9 GeV and = 60 GeV - the ones that best 
fit Coma astrophysical data - are too light to accommodate 
the PAMELA positron fraction (Adriani et al 2009) or the 
electron plus positron fiux seen by Fermi-LAT (Abdo et al. 
2009) and HESS (Aharonian et al. 2008, 2009b). 

In conclusion, the possibility of interpreting the origin of 
non-thermal phenomena in galaxy clusters with DM anni- 
hilation scenarios would require a low or intermediate neu- 
tralino mass and a cored DM density profile. If we then 
consider the multimessenger constraints to the neutralino 
annihilation cross-section, it turns out that such scenarios 
would also be excluded unless we introduce a substantial 
boost factor due to the presence of DM substructures. 
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Fig. 2. The radio halo spectrum of Coma and the besl^fits of 
three DM models with = 9 GeV (upper panel) for bb (solid 
with aV = 1.7-10"^'^ cm^ s"^ and Bo = 20 /xG) and (dashed 
with aV = 5.6 ■ 10"^^ cm^ s~^ and Bo = 2.2 fiG), = 60 
GeV (mid panel) for bb (solid with aV = 8.5 • 10"^^ cm^ s~^ 
and Bo = 0.5 ^G) and (dashed with aV = 9.0 ■ 10"^^ cm^ 
s"^ and Bo = 0.1 fiG), and with = 500 GeV (lower panel) 
with bb (solid with aV = 2.0-10"^* cm^ s"^ and Bo = 0.01 ^G), 
(dashed with aV = 1.9 ■ 10"" cm^ s"^ and Bo = 0.002 /xG) 
and (dotted with aV = 1.3 ■ 10"^*^ cm^ s"^ and Bo = 0.001 
/iG). Data from Thierbach et al. (2003). The results are obtained 
using B(r) oc nth{r) and NFW smooth DM profiles, without 
considering the effect of substructures. 
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Fig. 3. The radio halo spectrum of Coma and the predictions 
of DM models calculated by assuming Bo = 4.7 /iG and B{r) oc 
nth{'r)'^'^\ black lines refer to = 9 GeV with hh (solid line; 
aV = 6.0 X 10"^'^ cm^ s'^) and r± (dashed line; aV = 1.4 x 
10"^'^ cm^ s~^); red lines refer to = 60 GeV with hh (solid 
line; aV = 1.0 x 10"^"* cm^ s"^) and r± (dashed line; gV = 
9.0 X lO"^'^ cm^ s-i); blue lines refer to = 500 GeV with hh 
(solid line; aV = 7.5 x 10"^* cm^ s"^), r± (dashed line; gV = 
3.0x10"^^ cm^ s"^) and VF* (dotted line; aV = 1.1x10"^^ cm^ 
s~^). The results are obtained using NFW smooth DM profiles, 
without considering the effect of substructures. 
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Fig. 4. The SED of the DM-induced e.m. signals in Goma as 
predicted by the three DM models. Upper panel: = 9 GeV 
(solid: hh; dashed: t*); mid panel: My, = 60 GeV (solid: hh; 
dashed: t^); lower panel: My, = 500 GeV (solid: 66; dashed: r*; 
dotted: W ) The values of aV and Bo for the different models 
are the same used in Fig[21 The hatched regions in the X-ray 
domain represent data from PDS/BeppoSAX (Fusco-Femiano 
et al. 2004), HEXTE/RXTE (Rephaeh et al. 1999) and EUVE 
(Lieu et al. 1999). We also report the EGRET (Reimer et al 
2003), Fermi (Mori 2009) and HESS (Aharonian et al. 2009a) 
upper limits on Coma. The blue line is a thermal bremsstrahlung 
model for Coma with kT = 8.2 keV (Briel et al. 1992). The 
results are obtained using B{r) oc nth{r) and NFW smooth DM 
profiles, without considering the effect of substructures. 
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Fig. 5. The heating rate induced by DM secondary particles 
in Coma as predicted by the three DM models with = 9 
GeV, = 60 GeV and = 500 GeV for bh (solid), r± 
(dashed), (dot-dashed) compositions. We consider an NFW 
DM profile (green, blue and red peaking curves for 9, 60 and 500 
GeV respectively) and a cored DM profile (green , blue and red 
flattening curves for 9, 60 and 500 GeV respectively). The values 
of aV and Bq for the different models are the same used in Fig[2] 
The solid black curve shows the bremsstrahlung cooling rate of 
the intra-cluster gas at a temperature of 8.2 keV. The results are 
obtained using B{r) oc Uthir) and smooth DM profiles, without 
considering the effect of substructures. 
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Fig. 6. We show the constraints set by the radio halo spectrum 
of Coma on the annihilation cross section for the DM models 
considered in this paper. Upper, middle and lower panel refer to 
the 9, 60 and 500 GeV neutralino mass, respectively. We consider 
the case of an NFW DM profile with bb (solid red curve) and 
(dashed red curve), and the case for a cored DM profile with bb 
(solid blue curve) and (dashed blue curve). The lower panel 
contains also the cases of W* composition with NFW (dotted 
red curve) and cored DM profile (dotted blue curve). The upper 
curve bundles refer to a smooth DM profile while the lower curve 
bundles refer to the case in which substructures are considered 
as in Fig[l] The shaded bands in this figure represent the upper 
limit to the annihilation cross-section for the DM models derived 
combining the limits coming from a multi-messenger analysis 
(see text for details). The results are obtained using B{r) oc 
nth{r). 



